Two-dimensional 5 = 1/2 antiferromagnet on a plaquette lattice 



S.V. MeshkovB 

Laboratoire des Verres 
Universite de Montpellier II, CNRS UMR 5587 
Place E.Bataillon, 34095 Montpellier, France 

D. FoersterEI 

Centre de Physique Theorique et de Modelisation de Bordeaux 
Universite de Bordeaux I, CNRS URA 1537 
rue de Solarium, 33175 Gradiqnan, France 



On 
On 

u 

or 

< 

CO 

(N 



O 
o 



We consider a simplified model of the magnetic structure of the two-dimensional compound 
CaViOs in terms of interacting square plaquettes of spins with two distinct antiferromagnetic ex- 
change constants. We analyze the competition between two types of singlet ground states and 
the Neel ordered one in terms of respectively, numerical cluster expansion and nonlinear spin wave 
theory. The resulting phase diagram agrees well with known Quantum Monte Carlo results and sug- 
gests a first order transition between ordered and singlet ground states as a function of the exchange 
constants. 

PACS numbers: 75.10.-b, 75.10.Jm 



The recent experimental observation of a spin gap in the layered S = 1/2 antiferromagnet CaV^Og |Q| has opened 
a new and interesting perspective in two dimensional magnetism. We discuss the simplest model of the undoped 
I/") ' structure that consists of a square lattice of elementary squares or " plaquettes" , which we will refer to as C AVO 
lattice. The magnetic exchange energy within the plaquettes (Jo) and between the plaquettes (Ji) is given by 
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Here ij represent nearest neighbors on edges of a plaquette and between adjacent plaquettes, respectively. Addi- 
tional (frustrating) couplings are ignored, although they are believed to be necessary for quantitative agreement with 
experiments. 

I , The nature of the ground state is easy to understand in two limits In the limit of J\ -C Jo the plaquettes 
form resonating valence bond type singlet states, with an energy of —\J$ per plaquette, and weak bonds Ji serve 
as a perturbation. In the opposite limit of Jo <C Ji the interplaquette connections form singlets of energy — | Ji per 
dimer that are weakly interacting via plaquettes. This construction is qualitatively symmetric, but the plaquettes are 
somewhat "stronger", so that the critical "equilibrium" region is centered at Ji ~ |Jo- In this region, in addition to 
these two quantum singlet phases, the antiferromagnetically ordered Neel state could also be competitive. 

A variety of theoretical methods was applied to study the ground states of this model as a function of Ji/Jo but 
the results remain contradictory. Our purpose is to compare the energies of three candidates for the ground state 
of the model ([!]) and estimate the regions of their stability in terms of the ratio J\/Jq. Unlike other approaches 
that attempted to treat the system within a unified framework for all Ji/Jo, we choose the most quantitatively 
reliable approach for each individual phase. For the two singlet states we develop a numerical perturbation expansion 
(beyond the second order Q) in the coupling ratio Ji/Jo or Jo/Ji, whichever is smaller. The energy of the Neel 
state is estimated via the nonlinear spin wave approximation, which gives a lower ground state energy than the linear 
approximation reported in Q| . We also compare our results with ground state energies obtained by direct numerical 
diagonalisation of the Hamiltonian for different finite lattices of up to 24 spins and with frustrating/nonfrustrating 
boundary conditions. 



I. CLUSTER EXPANSION FOR THE GROUND STATE 



The idea of the cluster expansion is quite general and works for any model with finite range interactions provided 
the unperturbed Hamiltonian is a sum commuting blocks. For the model under consideration and in both limits 
Ji <C Jo and Jo <C Ji the perturbation is a sum of exchange interactions between nearest neighbors, and the zero 
approximation (the first or second term in Eq.([j]) respectively) is a sum over independent plaquettes or dimers. Thus 
in any given order n of the perturbation parameter (the smaller of the coupling ratios) the total correction to the 
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ground state energy of an arbitrary lattice cluster c can be reorganized into a sum over connected graphs that mark 
the interactions that were used and the unperturbed blocks that were touched: 

=Y,N c , g e g {n) ■ (2) 
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Here s g (n) stands for the contribution of the graph g in the order n and N c , g for the number of ways the graph g 
can be embedded in the cluster c. The same relation holds for the entire infinite lattice with properly normalized 
embedding numbers N x (g). 

The detailed analysis is rather tedious (see e.g. similar analysis || and references therein), and we only list the 
resulting sequences of graphs and embedding numbers. To be specific, consider the J\ expansion about the plaquette 
state, with four graphs contributing up to the 5th order 
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We draw the graphs as being embedded in the CAVO lattice, marking the Jo connections by bold lines. The Roman 
superscript of a graph denotes the lowest order of perturbation in which it contributes to Eq. (Q) . A connected graph 
contributes to order n if n interactions can be placed on the links of the graph (one or more per link) such that (i) 
each block is touched at least twice (otherwise the block cannot return to its singlet ground state), (ii) no two parts 
of the graph are connected by only one interaction (it would vanish by spin symmetry) . Once the contributions s g (n) 

(n) 

of all necessary graphs are known, the n-th order corrections Ec to the ground state energy of an arbitrary lattice 
cluster c can be calculated using the embedding numbers N Ct9 . One can, however, invert the procedure and recover 
the contributions of the graphs by solving the system of linear equations (0) for e g (n) in terms of the n-th order 

(n) 

corrections E c of an appropriately selected set of small clusters c. 

The most economic or "optimal" set of clusters is obtained when each cluster embeds some graph in the list of 
graphs exactly once. Note that the graphs are topological entities, so that the choice of optimal clusters is generally 
not unique. Once a choice has been made, there is a one-to-one correspondence between the graphs contributing to 
a given order and the finite clusters carrying the information on their contributions. Therefore we can use the same 
pictures and labels for optimal clusters as for graphs. For the plaquette expansion the embedding numbers N c ^ g of all 
graphs of (||) in all clusters of (||) are 

/ 1 
jupla _ 2100 
c -9 2 1 
\ 4 4 1 

with rows of the matrix corresponding to clusters and columns corresponding to graphs. 

The coefficients Ec are easily extracted from a polynomial fit of high precision ground state energies of the clusters 
in the list at several values of the coupling ratio. Having solved for e g (n) from Eq.(j|) we calculate (to the same 5th 
order) the ground state energy per site of an infinite lattice with embedding constants N^ a g = j (2 2 41) and obtain 
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So we find 5 orders of the perturbation series for the infinite lattice by diagonalizing only 4 finite and relatively small 
clusters of up to 16 spins by use of conventional Lanczos algorithms. 

In the case of the dimer expansion in the small parameter Jo/ Ji the unperturbed blocks are smaller (2 sites), so 
that we can reach 7th order of perturbation with clusters not exceeding 12 sites. The list of graphs/clusters contains 
13 entries 

l(")| 2( /y 'LL 3 (/y) _H 4( /y )[Zl 5 (/y) _tT" & {VI) \H 7( 1// )_hj- 
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in the same notations as Eq.([|), with the following matrix of embedding numbers 
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The embedding constants for the infinite lattice are N^™ — j (4 8 4 1 1 16 4 16 8 8 8 44), and the ground state energy 
per site in the dimer expansion is 
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II. NONLINEAR SPIN WAVES 



We apply the conventional Holstcin-Primakoff formalism, that parametrizes spins in terms of harmonic oscillators. 
This approach provides results of a remarkable precision for the s — 1/2 Heisenberg antiferromagnet on a square 
lattice (see || for review). The CAVO lattice is treated as a square lattice with nodes containing 4 spins. In the 
approximation of noninteracting spin waves, we find 
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r(Jo, Ji) = - (j + s(s + 1) - | J Tr{^l-7 2 (P, Jo, Ji)} 
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where the integral in p = (pi,p2) is over the square 2ir x 2tt Brillouin zone. This is analogous to the result for the 
simple square lattice except that 7(p, Jo, Ji) is now a 4 x 4 matrix with two exchange parameters 
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The next correction of order O(s ) is the average of quartic terms in the boson Hamiltonian in the linearly recon- 
structed ground state. Its direct evaluation in the Holstein-Primakoff formalism is rather tedious (see fl), but there 
exists a way to simplify the calculation. For bipartite lattices with two equivalent sub-lattices one can prove 

(S a S b ) = (s + Cabf 

for any pair of neighbor spins with c ao being a constant of order sr. In other words, every (S Q Sf,) is a full square in 
the nonlinear spin wave approximation. Therefore, the nonlinear spin wave result may be found by (i) separating the 
contributions from different types of nearest neighbors, (ii) completing the square for each of them and (iii) summing 
up the results. Thus for the model (0) we arrive at 
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where 7«(p) denotes 70 (p) 
substituted. 



7(p, 1,0) or 71 (p) = 7(p, 0,1). For the model under consideration s = 1/2 should be 
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III. DISCUSSION 



The results of our estimates for the infinite CAVO lattice are collected together in Fig.l. In these figures we used 
the energy of the classical Neel state ^classical (Jo, Ji) = — (2 Jo + Ji) /8 as a unit and we plot the rescaled ground 
state energy E = E(Jq, Ji) /^classical (Joi Ji) as a function of the reduced coupling J = J\/{Jq + Ji) using notations 
that are similar to those adopted in Q). 

For each of the two perturbative expansion we plot the energy in all computed orders starting from the 2nd. A 
tentative Pade extrapolation indicates closest singularities at respectively Ji/Jo ~ 1.4 for the plaquette expansion 
and J\/Jq ~ 1.05 for the dimer one, but the orders we reached are not sufficient to reveal the analytical structure 
of the expansions. However, both expansions appear to converge well in their dominant region, and we believe that 
the highest order of perturbation provides a good estimate for the ground state energy. Another fact to be noted is 
that the ground state energy per site £ rv b( J , Jo) = — 0.5499Jo following from Eq.(^) at the symmetricpoint J = Ji 
agrees perfectly with variatonal Monte Carlo calculations of the singlet state that give -0.5510 • Jo @ -0.5495 • J 

The precision of the 1/s expansion is generally less controlled than that of ordinary perturbation theories, as the 
small parameter is not obvious (see || for discussion). However, we see that at the symmetric point ( Ji = Jo) the first 
correction of order s (linear spin waves) gives ~ 45% of the classical Neel energy and the second, nonlinear correction, 
is ~ 10% of the first one. Such a convergence is only a little worse than in the case of the square lattice where these 
ratios are 32% and 8%, respectively, and we expect negligible further corrections. 

From Fig.l we see that each of the three states minimizes the energy in some region of Ji/Jo- Namely, the Neel 
state is stable in the interval 0.90 < Ji/Jo < 1.6, and the plaquette and dimer singlets are stable correspondingly 
below and above this region. This is in a perfect quantitative agreement with the Monte Carlo simulations ||, but 
not with the approximate treatments of [|l0| and jnj, where the Neel interval is considerably overestimated. Neither 
of our three curves for the ground state energy shows an anomaly in a reasonable vicinity of the intersection point. 
Thus we conjecture the occurrence of discontinuous (first order) transitions as a function of J\/ Jq, as a result of a 
direct competition in energy. A somewhat similar transition has been reported for a square lattice with additional 
frustrating couplings [ fl2|| . In fact, presently available results obtained by finite temperature simulations or mean- field 
type approximations cannot exclude this scenario. It assumes that the gap does not vanish on the border of the singlet 
regions, in agreement with recent perturbative estimates |1| extrapolated to the point Ji/Jo = 0.9. 

On the other hand, the fact that we use entirely different approaches for each region does not allow us to insist 
on the first order character of the transitions. The situation might be sensitive to minor quantitative changes due to 
higher order correction, so that the intersecting curves could turn out to be tangent. 

We found it interesting to complete the picture by direct Lanczos type diagonalization of finite clusters. The results 
for four clusters and their configurations are presented in Fig. 2. Three clusters are cut out from the infinite lattice 
with periodic boundary conditions imposed. The cluster of 4 plaquettes is nonfrustrated, whereas those of 5 and 
6 plaquettes are both frustrated. We have also considered another nonfrustrated cluster of 6 plaquettes designed 
"artificially" as an octahedron with vertices decorated by plaquettes. The energies are rather scattered due to small 
cluster sizes, but we observe that frustrated clusters tend to follow the singlet perturbation curves, while non frustrated 
ones follow the energy of the Neel state. To verify our cluster expansion we have computed the perturbation expansion 
for the 6 plaquette octahedron (to 5th order in the plaquette phase and to 6th order in the dimer one). The results 
shown in Fig. 2 suggest that the perturbation series are very precise in the regions of interest. 

The energy spectrum of the frustrated 5 plaquette cluster is rather unusual. Due to a special symmetry, two singlet 
states belonging to different representations intersect near J\ « 1.3 • Jo, while the lowest triplet state is continuous in 
Ji/Jq. Although this curious example may be quite special, we consider it as confirming the possibility of first order 
transitions on the infinite CAVO lattice. 
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FIG. 1. The rescaled ground state energy per spin E versus the reduced coupling J. Thin and solid lines represent, 
respectively, the linear and nonlinear spin wave approximations ([?]) and (^). Dashed ascending and descending lines represent, 
respectively, the perturbation expansions for the plaquette (Q) and dimer singlets starting from the 2nd order. The highest 
orders (5th for plaquettes and 7th for dimers) are bold-dashed. The two (hardly distinguishable) diamonds are variational 
Monte Carlo results M and Q. 
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FIG. 2. Energies of different finite clusters cut out from the CAVO lattice and that of an octahedron made up of plaquettes. 
Bold solid and dashed lines show the nonlinear spin wave and perturbative results on the infinite lattice. Dotted lines represent 
the perturbation expansion for the 24-site octahedron in the plaquette and dimer phase up to, respectively, 5th and 6th order. 
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